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Abstract 

The aim of this paper is to establish estimates of the lowest eigen- 
value of the Neumann realization of (iV + BA) 2 on an open bounded 
subset of M 2 Q with smooth boundary as B tends to infinity. We intro- 
duce a "magnetic" curvature mixing the curvature of dQ and the nor- 
mal derivative of the magnetic field and obtain an estimate analogous 
with the one of constant case. Actually, we give a precise estimate of 
the lowest eigenvalue in the case where the restriction of magnetic field 
to the boundary admits a unique minimum which is non degenerate. 
We also give an estimate of the third critical field in Ginzburg-Landau 
theory in the variable magnetic field case. 

1 Introduction and statement of main results 

Let fl be an open bounded subset of R 2 with smooth boundary and A G C°°(Q, 
We let : 

p = V x A 

and for B > and u e H l (Q) : 

4,nH= / \(iV + BA)u\ 2 dx 
Jn 
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and we consider the associated selfadjoint operator, i.e the Neumann real- 
ization of (iV + BA) 2 on fl. We denote by Xi(BA) the lowest eigenvalue of 
this operator. By the minimax principle, we have : 

\ 1 (BA)= inf Q ^ n{u) 



ueH 1 (Q) \\u\\ 2 

We first recall some properties of the harmonic oscillator on a half axis (see 
|DH93LIHM01] ). 

Harmonic oscillator on a half axis 

For ( G K, we consider the Neumann realization t) N ^ in L 2 (IR + ) associated 
with the operator 

(1.1) _^ + (t + e) 2 , V(t) N *) = {u G B 2 (R + ) : u'(0) = 0}. 

One knows that it has compact resolvent and its lowest eigenvalue is denoted 
/x(£) ; the associated L 2 -normalized and positive eigenstate is denoted by 
Uj: = «(•,£) and is in the Schwartz class. The function £ i— > admits a 
unique minimum in £ = £o an d we let : 

(1.2) e = M£o), 

uj (0) 

(1.3) d = -^i. 

Let us also recall identities established by |BS981 p. 1283-1284]. For k G N*, 
we denote by M k : 

M k = [ (t + ^ ) k \u^t)\ 2 dt. 
Jt>o 

(1.4) 

M = l, M 1 = 0, M 2 = ^, M 3 = y and = 3C iy /e~ . 

Let us state a result in the case where (3 is constant : 
Theorem 1.1 Assuming that (3 = 1, we have the estimate : 
X 1 (BA) = & B - C x n max 4B + O^ 3 ), 

where 

K max = max{A;(s), s G dQ} 

and k(s) denotes the curvature of the boundary at the point s. Moreover, the 
grounstate decays exponentially away from the points of maximal curvature. 
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Remark 1.2. 

This result was first announced by a formal analysis in [BS98J and rigorously 
proved in the case of the disk (see |BPT98| ). Let us also mention that in 
|LP99a| . an estimate at the first order was rigorously proved (see also |LP00| 



for the problem in IR 2 and For higher order expansion in the case of 

constant magnetic field, one can finally mention [dPFSOOl IHMOll IFH06al 
iFHflS] . 



Our aim is to obtain a similar result when the the magnetic field is not 
constant. We will assume that ft > on Q. We introduce : 

(1.5) b = in£ft and b' = iuift, 

U dVL 

and we assume : 

(1.6) Q b' < b. 
Estimate for the variable magnetic field 

Let us state a first (rough) estimate concerning the first eigenvalue : 

Theorem 1.3 Assuming that ft\gn admits a unique and non degenerate min- 
imum, we have : 

Xi(BA) = Q b'B + 0(B 1/2 ). 

Remark 1.4. 

The first term was obtained by many authors (cf. |LP99al IHM01| ) with a 
worse remainder estimate. Our assumption of non-degeneracy permits to 
find the optimal remainder 0(B 1 ^ 2 ) (the improvement occurs for the lower 
bound) which is crucial to establish tangential Agmon estimates (see Section 
4). 

■ 

Let us also state a tangential localization result of the first eigenfunctions : 

Proposition 1.5 (Tangential Agmon's estimates for u B ) Letu B be an 
eigenfunction associated with the lowest eigenvalue of the Neumann realiza- 
tion of (iV + BA) 2 . We have the control : 

J exp(aa(t(x))d(s(x))£ 1/2 ){|u B | 2 + B _1 |(iV + BA)u B \ 2 }dx < C\\u B \\ 2 , 

where x is a smooth cutoff function in a neighborhood of the boundary, t(x) = 
d(x,dQ), s(x) the curvilinear coordinate on the boundary and where d is the 
Agmon distance to the minimum of ft defined in Section 4- 
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Remark 1.6. 

This estimate improves the localization found in |HM01| by specifying the 
behaviour of ub near the minimum of j3. In Section 4 we also get tangential 
Agmon estimates for D s ub- All these localizations properties are essential 
to obtain the second correction term of Theorem 11.31 

■ 

Theorem 1.7 Assuming that (3\qq admits a unique and non degenerate min- 
imum in Xq, we have : 

X^BA) = & b'B + e 1/2 6 ,1/2 5 1/2 + 0(5 2/5 ), 

where 

e 1/2 = e 1/2 (x„) = -c(x„)c 1+ (| - e 0?0 ) if („)+eS" (^f^(-o)) 1 ' 

Remark 1.8. 

1. When (3\gn admits a finite set M. of non degenerate minima, we have 
the same expansion by replacing 61/2 by min9i/ 2 (x). 

2. Without assuming the non degeneracy of the minima, we believe that 
the conclusion of Theorem II .71 is true by replacing Q 1 / 2 by min Q 1 / 2 (x). 

3. The optimal remainder is certainly 0(B l / A ) as suggested by the upper 
bound. 

4. The computations for the upper bound suggest the following expansion 
of the n-th eigenvalue : 

\ n {BA) = e b'B + Q n 1/2 B 1/2 + 0{B 1/A ). 

where : 

e ?/2 = - K (x o) c 1+ (! - e ofo ) If^+pn-Der (^^)) 

5. In the variable case, the localization due to the curvature doesn't play 
a role anymore ; the effect of the curvature is small compared to the 
variation of the magnetic field. 

6. This expansion with two terms of the first eigenvalue could be general- 
ized at any order under the previous assumptions (unique and non de- 
generate minimum of f3\on) by using a Grushin approach (see |FH0 6aJ). 



4 



7. The case where the magnetic field (non degenerately) vanishes in O 
was treated in |KP02| . Moreover, the case where it non degenerately 
vanishes on the boundary remains open and should be an interesting 
problem. 

8. Theorems II .31 and II .71 are also sensible under the hypothesis of regular- 
ity of the domain. When the domain has corners (see |Bon05l Theorem 
1.2]) and with a variable magnetic field, the ground state is not nec- 
essarily localized near the points of the boundary where the magnetic 
field is minimum. 

9. The asymptotic behaviour in Theorems 11.31 and 11.71 is strongly depen- 
dent on the Neumann boundary condition we impose, as one can see 
in |Kac06l lKac07al rKac07b| . In particular, in certain cases, the local- 
ization is no more determined by the minimal points of (5. 



Constant magnetic field on the boundary 

In |Ara07l IAra06| , the case of the constant magnetic field on the boundary is 
treated. Nevertheless, this case is studied under a non degeneracy condition : 
it is assumed that the curvature of the boundary k admits a unique maximum 

at x = Xq and that the normal derivative admits a unique minimum at 

ot 

d(3 

x = Xq ; moreover, the minimum of — b'rz has to be non degenerate. Here, 

at 

we improve his result by using more generic assumptions ; in particular, we 
will see that the quantity to maximize is the "magnetic curvature" defined 
by: 

k(x) = C lK (x) + (W, - ^ ) l^{x). 
More precisely, our result is the following : 

Theorem 1.9 (Upper bound : constant magnetic field on dVL) When 
the magnetic field is constant on the boundary, we have the upper bound : 

\\BA) < e &'B-max jd/<*) " (y " ©o£o) ^(x^^B^+OiB 1 ^), 

where n{x)denotes the curvature of the boundary at x. 
Remark 1.10. 

1. The corresponding lower bound could certainly be obtained by the 
techniques of |FH08| . 
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2. Assuming the existence of a unique and non degenerate maximum of the 
magnetic curvature k, one could surely give an asymptotics at any order 
of X l (BA) and localization properties as for the constant magnetic field 
case (see |FH06a| ) which would improve the hypothesis of Aramaki. 



Organization of the paper 

In Section 2 and 3, we will prove the Theorem 11.31 and give the upper bound 
of Theorem 1 1 . 71 and of Theorem II .91 Then, we will see, in Section 4, that this 
first rough estimate gives information on the localization of the groundstates 
on the boundary near the mimimum of the magnetic field. In Section 5, we 
prove the lower bound of Theorem II . 71 thanks to a reduction to a degenerate 
case studied by S. Fournais and B. Helffer. Finally, we apply the previous 
results to give an estimate of the third critical field in Ginzburg-Landau 
theory. 

2 A rough lower bound 

In order to get the lower bound in Theorem 11.31 we use a localization tech- 
nique permitting the reduction to easier models. 

2.1 Partition of unity 

For each 0<p<~,-B>0,e>0 and C > 0, we consider a partition of 
unity (cf. |HM04p for which there exists C = C(Q, (3, e, C ) > such that : 

(2.7) J>f| 2 = lonfi; 

3 

(2.8) X)|Vxf| a <C7S^onn. 

3 

Each xf is a C°°-cutoff function supported in Dj n D,. Moreover, we may 
assume that there exists a ball Dj = Dj min whose center is the minimum of (3 
on the boundary and C B~ P for radius. We may also assume that the balls 
which intersect the boundary have their centers on the boundary and that 
those one admit eB~ p for radius. The radius of all the other balls is assumed 
to be B~ p . We will choose p, e and Co later for optimizing the error. We will 
use the following localization IMS formula (cf. (CFKS 86J) : 
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Lemma 2.1 

(2.9) q B A(u) = J2<lBA(xfu)-J2\\\ S7 xf\u\\ 2 , WueH\Q). 

j j 

So, in order to minimize qs a( w )j we are reduced to the minimization of 
1ba( v ), with v supported in some Dj. 

2.2 Estimates for the lower bound 

2.2.1 Study inside Q 

Let j such that Dj does not intersect the boundary. It is well known that : 

q B A(xfu)>B [ p(x)\ X fu\ 2 dx>bB [ \ X fu\ 2 dx. 
Jn Jn 

Having in mind (11.61) . these terms will not play a role in the computation of 
the asymptotics. 

2.2.2 Study at the boundary 

In the next paragraph, we introduce boundary coordinates. 
Boundary coordinates 

We choose a parametrization of the boundary : 

7 : R/(\dn\Z) -> dQ. 

Let u(s) be the unit vector normal to the boundary, pointing inward at the 
point 7(s). We choose the orientation of the parametrization 7 to be counter- 
clockwise, so 

det(Y(s),i/(s)) = l. 
The curvature k(s) at the point 7(3) is given in this parametrization by : 

7 "( s ) = k{a)v(s). 

The map $ defined by : 

$ : R/(|0fi|Z)x]O,i o [-> ft 
(s,t) 1 * 7(s) + tv(s), 

is clearly a diffeomorphism, when t is sufficently small, with image 
$(R/(|dfi|Z)x]0,to[) = e ft|cZ(x,dft) < t } = fV 
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We let : 

A x {s, t) = (1 - tk{s)) A($(s, t)) • y(s), A 2 (s, t) = A($(s, t)) ■ 

£( s ,t)=/WM)), 

and we get : 

8 i 3 -at^ = (i-tfc(aMM). 

Let j such that Bj intersect the boundary ; we have, with Vj = xf u an d 

Vj = Vj o $ : 

QBAivj) = J (I - tk{s))\{id t + BA 2 )vj\ 2 + (1 - tkis))- 1 ^ + BiO^fdsdt 

Approximation by a constant magnetic field on a domain with con- 
stant curvature 

Locally, we can choose a gauge such that 

A 1 (s,t)= [ (l-t'k(s))(3(s,t')dt', A 2 = 0. 
Jo 

We assume that the center of the ball Dj has the coordinates (sj, 0) and that 
the coordinates of the minimum are (0,0). We let : 

kj = k(sj), P(sj,0) = Pj and Akj(s) = k(s) — kj. 

We have : 
(2.10) 

(1 -tk(s))(3(s,t) = {l-tkj)Pj-t&kj(s)P(s,t) + {l-tk j )0{8,t)-Pj). 
We write : 

(2.11) A 1 (s,t) = A 1J (s,t)+R j (s,t), 
with 

(2.12) A ld (s,t) = (t-k~)P r 
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Control of the remainders 

Therefore, we are reduced to compare q B A with the quadratic form associ- 
ated with the Neumann problem on a domain with constant curvature (see 
[BPT981 lFH08l IFH06al IHMOlj ). For all A > 0, we get the inequality (with 
the Cauchy-Schwarz inequality): 

?ba(%) > (l-X)J(l- tk 3 )\d t v 3 \ 2 + (1 - tkj^Hids + BA^v^dsdt 

-C J Akjis^UdtVj] 2 + \(id a + BAJvjftdsdt 
B 2 f 

— — / \Rj(s,t)vj\ 2 dsdt. 

We apply the result of the constant magnetic field on a domain with con- 
stant curvature to get the existence of C > such that for all j such that 
Dj ndQ^tt (cf. |BPT98l Theorem 6.1]) : 
(2.13) 

J {l-tk^dtvtf+il-tk^^ids+BA^tfdsdt > (QofcB-dkjB^-C^VjW 2 . 



In order to control the remainders, we recall the Agmon estimates (cf. |Agm82 
IHMOll [FHflSl IFHOBap : 

Proposition 2.2 (Normal Agmon's estimates) Letu B be an eigenfunc- 
tion associated with the lowest eigenvalue of the Neumann realization of 
(iV + BA) 2 . We have the control the momenta of order n in the normal 
variable t : 

J t{x) n {\u B \ 2 + S _1 |(iV + BA)u B \ 2 }dx < C n B~^ \\u B \\ 2 . 

We choose p — \ (see Figure l2~T3l) and notice that |AjA;(s)| = 0(B^ 1 ^) 
(uniformly in j). So, there exists C > such that for all j : 



Akjis^ddtvjl 2 + \(id s + BA^v^dsdt 



< CBi\\vj 



r. II 2 



We let 



(2.14) a= 2^ ( °' 0) - 

Using the assumption of non degeneracy of the minimum, we can choose 
e > small enough such that 

(2.15) | S 2 < 0( 8 , 0) - 0(0, 0) < *-as 2 
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Figure 2.13: Partition of unity near the boundary 
for all \s\ < e . 

To estimate the other remainder, we will distinguish between three cases 

• \ s j\ > e o, 

• CqB~ 1 I a < \ Sj \ < e . 

Case 1 : j ' — j m i n 
As 

f(0,0)=0, 
we have, with (l2~Tnl) and fl2~TTD : 

\R^ n (s,t)\<C(t 2 + sH). 
Consequently, using Proposition 12.21 we get : 

/ \Rj - (s,t)i)j ■ \ 2 dsdt < CB'Hvj . II 2 . 

/ I Jmin V * / Jmin I — II Jratn II 

Taking A = B~ l l 2 , we deduce : 

q BA (v Jm J>(Q b'B-CB^ 



v, . || 2 . 
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Case 2 : \sj\> e 
We get : 



\Rj{8,t)\<C({8- 8j )t + 1*). 



Thus, we find : 



J \Rj(s,t)vj\ 2 dsdb < C(B~ 3/2 e 2 + B~ 



2 m\\ 2 - 



Moreover, there exists b" > b' such that for all \sj\ > e , we have : $j > b" . 
We take A = B^ 1 ^ 2 and deduce, using ( 12.131) and for B large enough, that 
for all j satisfying \sj\ > eo : 

QBx{vj) > Qob'BWvjW 2 . 

Case 3 : C Q B~ l l A < \ Sj \ < e 
We use the inequality : 



sup 

\s-s j \<eB- 1 / 4 ' 



<C sup \j3 -b'\ 



to find with fl2~TTD and ([2~T0D : 

f {Rjis^vtfdsdt < C(B~^ 2 e 2 sup \0 - b'\ + B- 2 )\\vj"' 



\s-Sj\<eB- 



As a consequence, we can write, with A = B l l 2 : 



q B A(vj) > (Qob'B + B{Q {/3( Sj ) - b') - Ce 2 sup \(3 - 6'|))||^| |2 

Is-s^eS" 1 / 4 

By non degeneracy, we have, for Co > 2e : 
(2.16) sup |/3-6'|<27 inf 

Indeed, we have, for all Co > 2e : 



a 



inf 1/5 - 6' I > - inf s 2 _ 

Is-^l^eB- 1 / 4 2 | s - Sj |<eB-l/4 



and 



3a 



3a 



sup \P-b'\<— sup s 2 < —( 8j + e^ 1 / 4 ) 2 . 



|s-Sil<efi~ 1/4 



{s-Sj^eB- 1 / 4 
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Thus, we get, for C > 2e : 



sup,,-,,,,,-., W -V\ (sj + eB^V ( 2eB~^ V < ^ 



infi,-,^-./. W ~ V\ ~ \>i-tB- 1/l J ' V sj-tB-'l* 
We deduce, for C > 2e : 

QBAivj) > {Qob'B + 5(00 - 27Ce 2 ) inf - &'|)||%|| 2 . 

We will further use that there exists c > such that for all C > 2e : 

QBAivj) > (e b'B + cB0( Sj ) - &')) ll^f • 
Indeed, we have, for all C > 2e : 

inf J/3-^>^fe)-&')- 
We find, for e > small enough : 

We conclude that : 

E^aK-)> (e ite - as 1 ' 3 ) £ N| a . 

j bnd j bnd 

Putting together this estimate and the estimate inside f2, we have the lower 
bound in Theorem II .31 

3 Models near a minimum of /3 and upper bounds 
3.1 Model operator 

We fix ko, ki and a > and we wish to study the quadratic form on the 

Hilbert space L 2 ((l — k t)dtds) defined, for u G C^(B ko ) by : 

(3.17) 

/k\ 
(l-tk )\d t u\ 2 +(l-tk y 1 \(-id s +Bt(l-—t+as 2 ))u\dtds J 



0<t< 



where B ko = M x 
associated operator is 



°> 2k 



(and by convention B = lx The self-adjoint 



^1 . , 2\\2 



;i - M) -1 9t(i - k t)d t + (i - tk ) 2 (-id s + 5t(i - yt + « s 2 )) 
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with Neumann condition on t = and Dirichlet condition on t = (if 
ko ^ 0). We first rescale the problem : 

t = B~ l /\ 

s = B-V**, 

and we are reduced to the operator on L 2 ((l — tkQB~ l ^ 2 )dtds) : 
(3.18) 

l-^VVl-^Wfl-^V7t-^t 2 + a t -^ 



B i/2 J 1 y B i/2 J 1 \ B 1 / 2 J \ 2B 1 / 2 B 1 / 2 B 1 / 4 

We make a change of gauge u i— > e^° Bl A(T u. Then, the operator defined in 

(13. 181) becomes : 

(3.19) 



+0O 

i=o 
13 



B i/2 J ^ y B i/2 J 1 \ B 1 / 2 J \ 2B 1 / 2 B 1 / 2 B 1 / 4 

3.2 Degenerate case : a = 

This case corresponds to the degeneracy of the minimum of the restriction 
of (5 to the boundary. In particular, we will prove Theorem 11.91 

3.2.1 Formal computation 

In order to have an upper bound, we first construct a formal quasimode. We 
make a Fourier transform in the variable s. Thus, we are reduced to the 
study of the family of operators on L 2 ((l — -^j^dt) : 

b*m* = -d - ^r'ftd - + (i - J&n* + « - ^rf^- 

We formally expand this operator in powers of B. 
Term in B° : 

H = -d 2 t +{t + if. 

Term in B~ l l 2 : 

H l = k d t -h(t + e)t 2 + 2k t(t + o 2 - 
We look for a quasimode expressed as : 



and a expansion of the first eigenvalue : 

+00 

j=0 

So, we have to solve 

H u = X u 

and, as we look for Ao minimal, we fix £ = £ 0) we deduce A = Go and we 
take Mo = U£ . Then, the next equation to solve is : 

H Q Ui + H-lUq = O Ml + AiUq. 

Thus, we deduce : 

(Ho - o )«i = (Ai - H^uo. 

To have solutions, the second member must be orthogonal to «o, so, using 
the formulas (11.41) . we get : 

Ai + ^y^Ck - e Uh - h) = 0, 

and we take : 

Ml = Ro(Xi - Hi)Uq. 

We let : 

Thus, ip is a good candidate to be a quasimode after truncation. 
3.2.2 Quasimode 

We write, in the initial coordinates (with b' = 1, for simplicity) : 

A 1 = A 1 + R, 

where 
with 

(3.20) = §^(0,0). 

Let us denote ip = u$ + i?" 1 / 2 ^! and notice that ^ is in the Schwartz class. 
As a quasimode, we take : 
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with x a smooth cutoff function supported in 



0> 2fc 



and p G]0, 4[ which 



will be choosen later to optimize the error. The Gaussian e s2b1/2 2p permits 
a localization near s = 0. We have : 

<?baM < J {I - tk )\d t u B \ 2 + (1 - tfo^K-ifl, + BA^UB^dsdt 

+C [ Ak{s)t{\d t u B \ 2 + \(-id s + BA x )u B \ 2 }dsdt. 



By noticing that there exists C > such that : 

\Ma,t)\ 

we get : 

Ak{s)t{\d t u B \ 2 + | (-10, + BA x )u B \ 2 }dsdt < CB l/4+f, \\u B \\ 2 . 

Let us prove the upper bound for the first term (the second can be treated 
in the same way). We have : 

^ MB = x '(t)^(5V 2 t)e- s2Bl/2 " 2p e^ Bl/2s +5 1 / 2 x(t)^(5V 2 t) e - s2Bl/2 - 2p e^ Bl/2 ^ 
Thus we get : 

|<W < 2\ X 'm(B l /H)\ 2 e- 2s2Bl/2 - 2p +2B\ X m\B 1 /H)\ 2 e- 2s2Bl/2 - 2p . 
Then, we find : 

J tAk(s)\d t u B \ 2 dtds < C j ts\x\t)\ 2 \i)(B l /H)\ 2 e- 2s2Bl/2 ~ 2p dtds 
+ CB f ts\x(t)\ 2 \^\B l / 2 t)\ 2 e- 2a2Bl/2 - 2p dtds. 



As if; is in the Schwartz class, we get : 

ts\x'{t)\ 2 \i,{B l /H)\ 2 e- 2s2Bl/2 - 2p dtds = 0(B~°° )\\u B \\ 2 . 



Then, we have after rescaling, for some C > independent of B : 

B j ts\x{t)\ 2 W{B l l 2 t)\ 2 e- 2s2Bl/2 ~ 2p dtds < CBB- l ' 2 B~ l /^\\u B \\ 2 = CB 1 ^ 
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Moreover, we have : 

J {I- tk )\d t u B \ 2 + (1 - tfco) _1 |H<9 s + BA 1 )u B \ 2 dsdt 
= J (1 - tfc )|<%w B | 2 + (1 - tko^K-ids + EZOmbI 2 ^ 
+3? |y (1 - tifco) -1 ^ 2 !^! 2 + 2B(-id s + BAi)u B .Ru B )dsdiJ . 
We get : 

feo^AiiW < (eoS + e^s^ + cs 1 ^!!^!! 2 , 

the crucial points being to estimate the term 

\d 2 s e- s2Bl/2 ~ 2 " \ e -° 2Bl/2 - 2p \ x (t) |X£ 1/2 t) \ 2 dtds 

by 0{B l / 2 - 2 P)\\u B f and the term $ {Bt + B 1 / 2 ^)d s (x{t)^{B l /H)e- s2Bl ' 2 - 2p ) 
by 0(S-°°)||w B || 2 thanks to the fact that M x = (cf. ([131)) and that ip is in 
the Schwartz class. Using that : 

\R{s,t)\ < C(t 3 + sH + st 2 ), 

we find : 



< CB^+quel 



3? |y"(l - tjfeo) _1 (5 2 |i2uB| 2 + 2B(-id s + B>4i)ub-Rmb)^J 
and finally with p = ^ : 

<Wn B ) < (e i? + e^B 1 ^ + ci^kii 2 . 

Thus, after replacing fci by its expression, the upper bound of Theorem 11.91 
is proved. 
Remark 3.1. 

It follows from the identities (11.41) that : 

^-0o£o = M 3 -£ o 3 >O, 

where M 3 = / (t + ^ofu^dt. This remark permits to understand how the 

Jt>o 

upper bound of Theorem 11.91 improves the one of Aramaki. 
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3.3 Non-degenerate case a > 
3.3.1 Formal computation 

We consider the operator H (cf. (13. 191) ) : 

Formally, we write : 

+oo 
3=1 

Let us look for a quasimode expressed as : 

+oo 

(3.21) U = J2 B ~ j/4U 3- 

3=1 

and a Taylor expansion of the lowest eigenvalue : 

+oo 

Af(5) = ^e j/4 ^/ 4 . 

3=1 

Here, we have : 

#o = ~d 2 t + (t + £o) 2 , 
H x = -2id s (t + Co), 

H 2 = k d t - d 2 s + 2(t + Zo)(asH - ^t 2 ) + 2k t(t + ^o) 2 - 
This leads us to solve : 

HqUq = XqUq. 

We write U as U = u (t)ipo(s) and, as we look for minimal, we take 
Ao = @o and uq > the associated normalized eigenvector. 
Then, we solve : 

^iZ7o + Wi = e ^i + Ai^o. 

We can take 0i/ 4 = by writing U\ = ui(t)ipi(s) with ipi = d s ipo and we 
find : 

(h - e )«i = 2i(t + e )«o- 

As Mi = (see (11.41) ). this last equation admits a unique solution U\ such 
that J t>0 UoU\dt = 0. 
Finally, we consider : 

h u 2 + h x u x + # 2 [/ = e t/ 2 + e 1/2 c/ . 
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Thus, we get : 

(H -e )U 2 = -H x U x -H 2 U Q + e y JJ Q = 2i{t + £ )u 1 d a ,/> 1 -H 2 U + e 1/2 U . 

Multiplying by u and integrating with respect to t, one applies the formulas 
f) 1 .41) and one solves : 

-(1 - 4/ 2 )<9 s Vo + ae s 2 Vo = (©1/2 + ^-^C x - (h - k )Q ^ V>o- 
where 

h= (t + £o)Ro((t + to)u )uodt. 
Jt>o 

This last integral can be rewritten by letting v = Ro((t + £o)^o) ; we have : 

(H - B )v = (t + Zo)u . 

By computing, we get : 

Idu 

-2^(^o) = - 
Using the identities of |FH08| . we find : 

i-4/ 2 = ^ = 3C lv / e;>o. 



After rescaling, we let : 

■00 (s) = e 2 v / 3^r 

and : 

(3.22) 1/2 = Q^ a = -^±h Cl + (h - k )0 Zo + v^ef yfa. 
3.3.2 Quasimode 

For simplicity, we assume b' = 1. We write : 

Ii = Ax + R, 

where 

M = t(l - ty + as 2 ) 
with a defined in fl2~T4l) and k x defined in fl3~20D . We let : 

u B (s,t) = X (t)U(B 1/4 s,B l / 2 t)e^ Bl/2s , 
18 



where U consists of the three first terms of (13.2 1 1) . We have : 

QbaM < J{1 - tk )\d t u B \ 2 + (1 - th)~ l \{-id s + BA^UB^dsdt 
+C [ Ak(s)t{\d t u B \ 2 + \(-id s + BA x )u B \ 2 }dsdt. 



Moreover, we have : 

J(l- tk )\d t u B \ 2 + (1 - tko^H-ids + BAJuBfdsdt 
= J (1 - tk )\d t u B \ 2 + (1 - tko^K-ids + BA^u^dsdt 

+3? |y (1 - tkQY l {B 2 \Ru B \ 2 + 2B(-id s + BA 1 )u B Ru B )dsdt^ . 

Using that U is in the Schwartz class, we get : 

QkoM^sM < (e B + Q^' a B^ 2 + C)\\u B f. 
Moreover, we have : 

l^i-^xl < C(sH + st 2 + t 3 ). 

So, we get : 

3? |y"(l - tk Q )- 1 (B 2 \Ru B \ 2 + 2B(-id s + BA 1 )u B Ru B )d8dtj < CB 1/4 \\u B 

Finally, we find : 

q B x{u B ) < (Q B + Q^B 1 ' 2 + CB l ' 4 )\\u B \\ 2 . 
In particular, we have proved the upper bound in Theorem 11.31 

4 Tangential Agmon's estimates 

We first observe that, for $ a real Lipschitzian function and if u is in the do- 
main of the Neumann realization of (i\7 + BA) 2 , then we have, by integration 
by parts : 

$t((iV+BA) 2 u,exp(2B 1/2 <$>)u) = q BA (exp(B 1/2 $)u)-B\\ | V$| exp(B 1/2 $)u\\ 

Taking u = u B an eigenfunction attached to the lowest eigenvalue X 1 (BA) , 

we get : 

(4.23) 

A^SAJIIexp^ 1 / 2 ^)^!! 2 = gBA (exp(5 1 / 2 $)u B ) - B\\ | V$| exv{B 1/2 $)u B \\ 2 . 
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4.1 Tangential Agmon's estimates for ub 

We now use the lower bound found in Section 2 ; more precisely, for all e > 0, 
there exists c > and C > such that, for all Co > sufficiently large, there 
exists C > s.t for all u in the form domain of qBA '■ 



q B A(u) > (bB-CB^J^WXj 



j int 

2 



+ (Oob'B + ciPis^-b'^WxjU 

j bnd,j^j min 

+(e b'B-C'B 1 / 2 )\\ Xj 



u\\ 2 . 



We choose u = exp(B 1 / 2 &)itB ; we recall that, by Theorem 11.31 we have the 
upper bound : 

X 1 (BA) < <d b'B + CB 1/2 . 
Using these estimates in (14.231) . we find the inequality by dividing by B : 

J (CB-W + \V^)\ Xjmm eMB l,2 ^)u B \ 2 > 
j (c0(sj) ~ b 1 ) - CB- 1 / 2 - \V^\ 2 )\x^MB 1/2 ^)u B \ 2 dsdt. 

j bnd 
3^3 min 

We choose 

$ = acid(s), 

where d is the Agmon distance associated with the metric (/?(s, 0) — b')ds 2 
i.e : 

M 

d{s) = / (P(a,0)-b') 1/2 da. 



o 



On Da . we notice that 

Jmvn " 

|V$| 2 < CB- 1/2 . 

Then, for j ^ j m in, we consider the quantity : 

c(P( Sj ) - b') - CB- 1 / 2 - al0(s) - b'). 

For e > and a.\ small enough, there exists d > such that for j such that 
\sj\ > eo and B large enough, we have : 

c0( Sj ) - b') - CB- 1 ' 2 - a\0{s) - b') > d. 
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For C > 2e, there exists c" > such that for j ^ j min and \sj\ < e and B 
large enough, we have : 

c0( Sj ) - b') - CB~ l/2 - a{{P{s) - b') > c"B~ l/2 . 



Indeed, due to the non degeneracy, we have (12.161) . Thus, we get C > and 
B > such that for all B > B : 

V f \ Xj expiB^usl 2 <C [ | eMB 1,2 $>)u B \ 2 . 

jb nd J J\s\<c B-m 

We deduce Proposition 11.51 and have the following corollary : 

Corollary 4.1 For all n £ N, there exists C > such that for all B large 
enough : 

[ s 2n {\u B \ 2 + B _1 |(iV + BA)u B \ 2 }dx < CB~ n ' 2 [ \u B \ 2 alx. 
Jn Jn 

4.2 Agmon's estimates for D s ub 

We consider a partition of unity as in (I2.7B . We have the formula (I2.9p and : 

q B A{u)>Y,<lBA{xfu)-CB l l 2 ^ 2 



it . 



We use (OHjl . We have 

A^B) < Q b'B + CB 1/2 . 
Thus, we get, using the inequalities of the previous section : 

Qba {x imin e B ^u B ) + ( Q ob'B + c((3( Sj ) - b')B) \ Xj e Bl '^u B 

+bBj2\x 1 e Bl/2 "'u B 



j int 



< (Q b'B + CB 1 ' 2 



where $ = a\d{s). 
We have the control : 



#1/2$ 

e u B 



B J \w^e Bl/2,s> u B \ 2 + CB^ 2 \\u B \ 



B 



„ Ti ^1/2$ 



< CB 1 ' 2 



n 
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< CB 1 / 2 / |u B | 2 dx, 
Jn 



and we deduce, for a.\ small enough : 

^A(Xi mi „e Bl/2 *n B ) - Q b'B 



Y, . e Un 



< / \U B 



We introduce : 
(4.24) q apP {v) 



[i-ht)\d t vf 



t>o,sen 



(1 - kot)- x \{Bt + 5as 2 t + B 1/2 £ -D s - B^-t 2 )v\ 2 dtds. 

If we write : 

A 1 (s,t) = J(l-t'k(s))p(s,t')dt', 

we have : 

(1 - tk(s))(3(s, t) = (1 - tfci) + as 2 + 0{t 2 + st + s 3 ), 

and thus : 
(4.25) 

Then, by the Cauchy-Schwarz inequality, we have for all A > : 

B 2 

Qba{v) > (1 - \)q apP {v) - — \\Rv\\ 2 . 

For instance, we can estimate §(st 2 ) 2 \v\ 2 . Using the tangential (cf. Proposi- 
tion 11.5ft and normal Agmon estimates and letting : 

v — Ya . e «r, 



Ai(s, t)=t - yt 2 + as 2 t + 0(t 3 + st 2 + sH). 



B 2 / s 2 t A \v\ 2 dsdt < CB 2 B- 1 l 2 B- 2 \\v\\ 2 . 



we have : 



In the same way, we control the other remainders and by choosing A correctly, 
we get : 



(4.26) 



Qba(v) >q apP {v)-CB l / A J 



Using the Cauchy-Schwarz inequality and again the Agmon estimates, we 
find : 

|2 



q apP (v) > (1 - B^ 2 ) q 2 app (v) - CB 1 ' 2 J 
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where 



7 2 ( 

iapp\ 



1 - k t)\d t v\ 2 + (1 - k ty l \{Bt + B 1/2 £o - D s )v\ 2 dx. 



t>0,s£l 



Making a Fourier transform in the variable s and letting w = v, we have 



[ (1 - k t)\d t w\ 2 + (1 - kot^KBt + B 1 ' 2 ^ - a)w\ 2 dtda. 
Thus, we get (see |FH081 Chapter 6, Prop 6.2.1] or [HM01L Section 11]) : 



q 2 app (v) > Qob'B J 



|2 _|_ ^1/2 f 1 " (£o) 



< C / |«b 



Consequenlty, we get the upper bound : 
We deduce the following proposition : 

Proposition 4.2 (Tangential Agmon's estimates for D s ub) With the pre- 
vious notations, there exists C > and a± > smc/i i/iai /or a// B large 
enough : 

2 



dx<CB 1/2 / I 



e a 1 B 1 ^ X (t(x))d(s(x)) DsUB 



where x is a smooth cutoff function supported in [— to, to]. 

Corollary 4.3 For all n G N, i/iere existe C > sttc/i i/mi /or ai/ B large 
enough, we have : 

[ X (t)s 2n \D s u B \ 2 dx < CB l/2 ~ n/2 [ \u B \ 2 dx. 
Jn in 

Remark 4.4. 

The tangential and normal Agmon estimates roughly say that \ub\ has the 
same behaviour as e~ as2Bl/2 uo(B 1 ^ 2 t). 



5 Refined lower bounds 

In this section, we prove the lower bound in Theorem 11.71 We consider a 
partition of unity as in (12.70 with p = - — i] for i] > 0. We have : 

9ba(«) > - CB^-^WuW 2 . 
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5.1 Control far from the minimum 

Let us first recall some the estimates we have proved. For j such that Dj 
does not intersect the boundary, we have : 



QBAiXju) > bB J \xju\ 2 dx. 



For j such that Dj intersect the boundary and j ^ j m in, we notice that, for 
B large enough : 

QBAiXjU) > Q b'B J \Xju\ 2 . 

5.2 Reduction to a model near the minimum 

Using the inequalities of the previous section, we get : 

QbaM > &ob'B J2 Wx^bW 2 + q B A(x jmin u B ) ~ C B l l 2 ^\\u B \\ 2 . 

By the normal and tangential Agmon estimates, we have proved in (14.261) . 
with ( 14.241) . ( 14.251) and the Cauchy-Schwarz inequality : 

lBA(Xj min UB) > q apP (Xj mm u B ) - CB l,A \\u B \\ 2 . 

In order to make the term in as 2 t disappear, we make the change of variables : 

t = A(s)r, 

where \(s) = (1 + as 2 )~ l l 2 ; we have 

dr _ _ dr 
(5.27) d s v = ^-d r v + d s v, d t v = -7r-d T v, 

os at 

where v denotes the function v in the variables (r, s) and we are reduced to 
the form : 



q apP (v) = J |(l - fc rA(s))|d 7 



2 



hr 2 



+ (l-fcorA(s)) 1 \{Bt + ^X(s)B 1 / 2 -X{s)D s -B^X{s) 



+ arsA(s) 3 D r )w| 2 |A(s)" 1 c/rc/s, 



where we have omitted the tilde. Noticing that s 2 = 0{B 2p 1 ^ 2 ), on the 
support of v = Xj m in u B, we make the approximations in L 2 : 

-X(s)D s v = -D s v + O (s 2 ) D s v, 
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T 2 X(s) 3 V = T 2 V + O (s 2 T 2 ) V, 

s\(s) 3 tD t v = stD t v + 0(s 3 r)D T v. 

We first find : 

q^p(v) > / {(1 - rk )\d T v\ 2 



+ (1 - rk )- l \ (Bt + £ \(s)B 1/2 - X(s)D s - B^\(s) 3 + asXisfrD^v^jXis^drds 

-C J &\{s)T{\dM 2 +\{Bt + ^\{s)B 1 ' 2 - \{s)D s - B^\{sY 
+ as\(sfrD T )v ] \'}\(sY l dTds, 
where 

AA(s) = A(s)-A(0). 
Let us consider the second term : 

J AX(s)r{\d T v\ 2 +\{Br + ^ X(s)B 1 / 2 - X(s)D s - B^X(s) 3 
+ asA(s) 3 rD r )t;| 2 }A(s)^ 1 rfrds. 
Coming back in the variables (t, s), this term becomes : 

/ ^f^ 1 ^' 2 + 1(5(1 + ^ + ~ ° s ~ 5 ^)^ 2 }^ s - 

Thus, the Agmon estimates give a control of the second term of order 0(1). 
Then, by the Cauchy-Schwarz inequality, the Agmon estimates (for ub and 
D s ub after having come back in the variables (s, t)) and using the same kind 
of analysis as in (14.261 ) , we have : 

J {{I- rk )\d T v\ 2 + (1 - rfco)- 1 ! (Br + ^ X(s)B 1 / 2 - X(s)D s - B^fx(s) 3 
+ asX^sfrD^v^jX^s^drds 

> J {(1 - tA; )|<9 t *;| 2 + (1 - Tk )- 1 \(Br + £ A(s)5 1/2 - D s - B^-)v\ 2 }X(s)- 1 dTds 
-CB l ^\\v\\ 2 . 
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We have finally, with v = Xj min u B '■ 
(5.28) 

Qba{ub) > Qob'B ^2 WxjUbW 2 



3 ¥^3 rain 

/7 2 
{(1 - rk )\d T v\ 2 + (1 - rko)- 1 ] (Br + ^ X(s)B 1 / 2 - D s - B-^v^X^drds 

- CB l ^\\u B \\ 2 - CB l ' 2 ^\\u B \\ 2 . 

Moreover, thanks to the exponential decrease of ub away from the boundary 
(normal Agmon estimates) , we can replace Xjmin by a smooth cutoff function 
such that 

suppXi_ C {0 < t < S-V^andlsl < E" 1 / 4 ^'}, 

that is we assume Xj min is supported in rectangles rather than balls ; the 
reason is technical and will appear in the next section. 

5.3 Lower bound for the model 

f s 

So, we are reduced, after the rescaling r = , s = , , to the study of : 

B I B I 

q m od(u) = [ {(1 - -^)\dM 2 
JT>o,sm k a 

+ d " Wnr\- + tMB-v>i) - ^ - Jt^H 2 }! 1 + 1^)" 2 **' 

Reduction to the euclidean measure 

In order to make disappear the measure (1 + Y^ 2 ■> we make the change 
of function defined by : 



-2 \ 1/4 

as N 



we have : 

^'O 7 Ma 1 2 



q m0 d(u)= [ {i 1 - -5m)\ d M 

+ (1 - + - ^ - - 5^H' 
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Term in s 

We want to make a Fourier transform in the variable s to be reduced to a 
problem on a half axis, but the term ^ X(B~ 1 ^ 4: s) is annoying ; that is why we 
make it disappear with a change of gauge. We write : A(£> _1 / 4 s) = 1 + t\b(s) 
and we make the change of gauge v i— > v = ve~^ 8 \ where 

to be reduced to : 

Jt>0,s€R k a 

+ d - ^r'K* +fo-^i- ^fjsi'Jdftu, 

where it = (Xj min Wfl)(S 1 ' 2 f, i? 1 / 4 ^). We make a Fourier transform in the 
variable s and we are reduced to a half axis problem in the normal variable : 

*<»> = jfjl - + d-|^)- I l(r + €c-^ 7 i- 5 ^ 7sf *»)»|»«*f, 

with w = v. 



Model on a half axis 

We can apply the same kind of analysis as in [FH081 Chapter 6, Prop 6.2.1] 



or in [HMOll Section 11] to get the lower bound ; there exists C > such 

that for all B large enough : 

(5.29) 

qn (w) > (e + (e** + ^. 2 )r'/ 2 - cir J m 2 (1 - 



Remark 5.1. 

In |FH08| . the fact that the magnetic field is constant permits to be reduced 

"'■o,k 
/2 



to the case fc = k\ — 1, thus Qy'^ 1 = 



Let us just recall the main ideas of the proof. We consider first the (formal) 
operator on L 2 ((l - -^)dr) : 

r / m t-i k r , d k f d _ k f . 2 , A ., a , f 2 
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Then, we formally expand this operator in powers of B and, for |er| < MB 11 , 
with rf > small enough : 

f)(a, B) = f) + B~ x ^fy v + £T 1/2 J) 2 + 0(£T 3 / 4+3?? ), 

where 

d 2 

^0 = "^2 + ( f + ^) 2 , 

^ = -2(f + e )a, 

^2 = - fcif 2 (f + &) + 2fc r(f + e ) 2 + a 2 . 

Thus, as in Section 13.3.11 we compute a quasimode and obtain for some ip : 
\M* 7 B) - (A + X.B-y + X.B-^m^^^ = 0(B-^). 

Finally, we can prove that the previous operator admits only one eigenvalue 
strictly less than 1 thanks to a comparison with the harmonic oscillator on 
a half axis and, applying the spectral theorem, we get the bottom of the 
spectrum given in (15 .291) (the values of a such that \a\ > MB V provide 
higher energies thanks to the non-degeneracy of £ h- > //(£) near £ )- 

Return in the initial variables 

Applying the Parseval formula, we get : 

q mo d(u) = q^d(v) > 

t >0 

+ B -i/2^M_ I ^2 f 1 _ rko\ dfd§ _ CB-^\\u\\ 2 . 

f >0 

We have : 

\Dgu\ 2 = \{D g -<f>'{a))v\ 2 . 
As \<f>'{s)\ < Cs 2 B~ l l 2 < CB' 1 ' 2 ^ on the support of v, we get : 

\D s v\ 2 > (1 - B-^)\D s v\ 2 - B~ 1/4+n \v\ 2 . 

Moreover, we have : 

OtS 

DsV = T^T72"^( s )~ 3m + fB{s)D 3 u. 
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We deduce : 

|£>^| 2 > \D s u\ 2 ~ CB-V^Qu] 2 + \D s u\ 2 ) - B-V 4+r >\D 
Recalling that 



,^l 2 . 



g \ 1/2 

dfds = ( 1 + -gq^ \ dids, 

s 1/2 

2 I i aS \ I |2 



where t = B~ 1 / 2 t, and with the tangential Agmon estimates, we get : 
HA^H 2 < C\\u\\ 2 , \\D s u\\ 2 < C\\u\\ 2 

and 

q m od(u) > Q \\u\\ 2 + Q^B-^WuW 2 
(5.30) + (| {ae |sn| 2 + ^M|Z^| 2 ds} (l - J^) IT 1 " 

where u(t, s) = u(f, s) and, thanks to the Agmon estimates, we have re- 
placed DgU by DgU and f by t by noticing that D§u = D§u + ff-DfW and 
A(s J B^ 1 / 4 )f = f. We recognize the quadratic form of the harmonic oscillator 
and we have : 



q;Bo|sm| H — as > > y ^ / M ds. 

We take V — ^ an d the lower bound of Theorem 11.71 follows from (15.301) . 
(15.281) and (11.41) after having noticed that the estimates of Agmon give : 



X jmin u B \ 2 dx= (l + 0(e ci * )) / M 2 dx. 

6 Estimate for the third critical field of the 
Ginzburg-Landau functional 

In this section, we give an estimate of the third critical field of the Ginzburg- 
Landau functional in the case where the applied magnetic field denoted by 
(3 admits a unique and non degenerate minimum on the boundary of VL. 
The constant magnetic field case has already been studied in details (see 
|FH06bl lLP99al ILP99H ILPOO] ). 
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Recall of properties of the functional 

The Ginzburg-Landau functional is defined by : 

W,A)= / {|(iV+aKA)^| 2 -K 2 |^| 2 +^^| 4 }^+W 2 / |VxA-/?| 2 rfx, 
Jn * Jn 

for ip G H\n,C) and A G fl^„(fi,R 3 ) where 

H^ iv (n,R 3 ) = {A e H\n,R 3 ) : div(A) = Oinfi, A ■ v = Oondtt}. 
We assume moreover that 

P = V x F. 
Then, we introduce the critical fields : 

Hc 3 (k) = infjcr > : (0, F) is the unique minimizer of G k ,ct}, 

Hc 3 (k) = inf{<7 > : (0, F) is the unique minimizer of Q K y for all a' > a}, 
H_c 3 ( K ) = m f{cr > : (0,F) is a minimizer of(? Kj(J } 

and 

7?cI{k) = sup{a > | A 1 (kctF) < k 2 }. 

We have 

H C3 (k) < H Cs (k) <H Cs (k) 

and 

H 1 ° c z {k)<HcM- 
We can prove the following result (cf. [FH08J) : 

Theorem 6.1 Let Vt be a bounded, simply connected domain with smooth 
boundary and suppose that the applied magnetic field (3 satisfies 

< Q b' < b. 

Then, there exists k > such that for all k > k : 

H C3 (k) = H 1 °1(k). 

Furthermore, if B i— > X 1 (BI) is strictly increasing for large B, then all the 
critical fields coincide for large k and are given by the unique solution H of 

\\kHF) = k 2 . 
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Estimate of H C3 {k) for large k 

Noticing that B i— > X 1 (BF) is strictly increasing for large B (it is due to the 
exponential decrease of the first eigenfunctions away from the boundary, still 
true in the case of variable magnetic field ; see [FH081 Chapter 9, Section 
6]), we deduce the following theorem : 

Theorem 6.2 Let Vt be a bounded, simply connected domain with smooth 
boundary and suppose that the applied magnetic field (3 has a unique and non 
degenerate minimum on dVL and that : 



< Q n b' < b. 



Then, we have : 



& /1/2 %! + o(«r 7 / 20 ). 



b'Q ef 
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